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Abstract 

The two-dimensional interpolating readout concept of the Virtual-Pixel detector 
(ViP detector) goes along with an enormous reduction of electronic channels com- 
pared to pure pixel devices. However, the special concept of the readout structure 
demands for adequate position reconstruction methods. Theoretical considerations 
reflecting the choice and the combination of linear algorithms with emphasis on 
correct position reconstruction and spatial resolution are presented. A subsequent 
two-dimensional coordinate transformation further improves the image response of 
the system. Measured images show how far the theoretical predictions of the simu- 
lations can be verified and to what extent they can be used to improve the position 
reconstruction. 
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1 Introduction 

In general, position sensitive pixel detectors feature the advantages of par- 
allel and asynchronous readout which means i.e. a proportional increase of 
the maximum incoming rate with the sensitive detection area. However, these 
pure pixel devices suffer from the enormous number of electronic channels 
leading either to high costs and effort or essential size limitation of the sensi- 
tive area. In contrast to that, interpolating readout concepts which can cover 
large areas with a position resolution comparable to pure pixel devices can 
also be carried out with parallel and asynchronous readout. They combine the 
advantages of pixel detectors with a reduction of electronics and costs at the 
same time. For example the interpolating readout can be realised by strips or 
pads [1,2,3,4]. We have developed a truly two-dimensional interpolating read- 
out concept based on resistive charge division [5]. Since this concept requires 
to count single photons the signal of each event has to be amplified, e.g. by 
gas amplification. For that purpose e.g. micro pattern devices like GEM ^ [6], 
micromegas ^ [7] and CAT ^ [8] can be applied. We have combined the readout 
structure with (optimised) MicroCAT gas gain structures [9,10] and a triple- 
GEM configuration [11]. In combination with the triple-GEM configuration 
this readout concept has already proven its possible application in biological 
and chemical diffraction measurements [12]. In principle this resistive interpo- 
lating readout structure can also be applied with vacuum devices like micro 
channel plates. 

Beside the advantages of a good time resolution [13] and high rate capability 
this two-dimensional interpolation concept demands on the other hand for 
a sophisticated choice of position reconstruction algorithms and correction 
functions. 



2 Detector setup and working principle 

Primary charge generated in the gas conversion region by incoming photons 
(typical energy range: 5 — 25 keV) via photoelectric effect is led by an homoge- 
neous electric drift field towards the gas gain region where it is multiplied in 
an avalanche process. This charge cluster (a so called event) hits the readout 
plane, designed as a two-dimensional resistive structure, where the sensitive 
area is subdivided into 7x7 square, 8 x 8mm^ sized cells (Fig. 1). Each cell 
is provided with two different surface resistivities in order to obtain both an 

^ GEM = gas electron multiplier 

^ micromegas = micro mesh gaseous structure 

^ CAT = compteur a trous 
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Fig. 1. Schematic of the two-dimensional interpolating resistive readout anode 
printed on a ceramic substrate. The gas gain structure (e.g. GEM or MicroCAT) is 
mounted above the anode. The structure is subdivided into 9x9 cells, whereby the 
inner 7x7 cells, corresponding to a sensitive area of 56 x 56 mm^, can be read out. 
A detailed magnification of one single cell of the interpolating readout structure 
on a ceramic substrate is illustrated on the right hand side. The charge, generated 
by a photon, is collected by readout nodes situated at the cell corners (represented 
by the black circles). The robust readout structure itself is absolutely insensitive 
against sparking or uncontrolled discharges. 

almost linear charge division behaviour in both coordinates and a good electri- 
cal shielding to avoid too large charge distributions over several cells. The high 
resistivity centres of the cells are surrounded by low resistivity narrow border 
strips (s. detail in Fig. 1). This leads to an almost perpendicular projection of 
the charges in x- and y-direction onto the low resistivity strips. The surface 
resistances of the high resistivity cell centres and the low resistivity strips are 
typically Rl = 100-1000 kf^/D and R2 = 1-10 kQ/D, respectively. Each cell 
is read out at its four corners by low input impedance readout nodes. The 
charges collected at the particular nodes are used to determine the event po- 
sition within the cell by means of a suitable algorithm. All 64 readout nodes 
are equipped with charge sensitive amplifiers and 12- bit FADCs which are 
linked to a dedicated DAQ system which transfers the preprocessed data for 
visualisation to a PC [14]. 



3 The choice of reconstruction algorithms 

The first subsection describes the simulation tool which provides the charge 
distribution within the readout structure. These theoretical data will be used 
in the following to study the partly unfavourable behaviour of different linear 
algorithms concerning position reconstruction and spatial resolution. After- 
wards, the treated algorithms are combined to one more complex algorithm 
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which unifies the advantages of the individual algorithms. 



A remark concerning the notation: The expressions for the derivatives dx/du 
and duX are used synonymous. 



3.1 The diffusion simulation model 



Due to the surface resistances and the underlying distributed capacitance the 
readout structure behaves like an integrating i?C-element resulting in a tem- 
poral broadening of the input signal. The simulation of this dynamic charge 
diffusion process on the readout structure has been realised by means of a 
numerical computation programme, taking 4 x 4-cells into account in order 
to include charge flow into adjacent cells (boundary effects). A very detailed 
derivation of the diffusion model and its numerical solution can be found in 
Ref. [15]. In the diffusion model a space and time dependent driving current 
can be impinged at every grid point. In the simulation, presented here, we 
assumed a typical two-dimensional gaussian-like transverse extension of the 
charge cloud of cr = 200 /xm (mainly caused by the transverse diffusion of the 
primary electrons in the conversion gap). 

The 4 X 4-cell model provides 161 x 161 grid points which correspond to a 
homogeneous grid point spacing of Ax = Ay = 200 //m. The low resistivity 

border strips have a width oi w = 200 ;um, hence they are represented by 
one single line of grid points. If not stated differently the readout nodes are 
represented by one grid point which corresponds to 200 x 200 /xm^. Since the 
readout nodes can be well approximated as ideal drains [15] the potential at 
the nodes is equal to zero at all times. 

This time-depending diffusion simulation is able to provide the currents In{t) 
and the charges Qn{t) at every particular node n for all times t. These simu- 
lated charges Qn{t) can now be used to reconstruct the event position (u, v) by 
means of a suitable algorithm. Depending on the algorithm the reconstructed 
position (m, v) will differ up to a certain extent from the (true) impact posi- 
tion (x, y). Since this simulation model has proven that its results are in good 
accordance to the measurements [15] it is used in the following to obtain the 
relationship between the two spaces (-u, v) and (x, y) to optimise the applied 
reconstruction algorithms. For that reason all further steps are directly based 
on the predictions of the diffusion simulation. 
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Fig. 2. This figure shows the node- and the cell-numeration. The origin of the 
coordinate system for the 6- and the 3-node algorithm is located at node 4. For the 
4- node algorithm the origin is chosen in the centre of cell 1. The cell size is denoted 
with g. 

3.2 The linear reconstruction algorithms 

Only in the theoretical case of R1/R2 = oo the charge division in two dimen- 
sions is exactly linear. The required high rate capability limits Rl upwards 
since large values of i?l lead to stronger temporal charge diffusion and hence to 
a slower charge collection at the readout nodes. On the other hand the paral- 
lel resistive noise which determines the spatial resolution significantly restricts 
R2 downwards. Therefore, only finite resistivity ratios of R1/R2 will be re- 
alised in practice. For that reason any linear reconstruction method produces 
approximate positions only. 

All our attempts introducing non-linear terms of the collected charges Q„ in 
a position reconstruction algorithm were not successful under the boundary 
condition of correct reconstruction at certain symmetry points of a cell. The 
discontinuity caused by the low resistivity cell borders which have been in- 
troduced to prevent the charge from leaving the cell and thus to improve the 
high rate behaviour obviously complicate this attempt enormously. Therefore 
these non-linear approaches are not further considered. 

The simplest and most obvious linear reconstruction methods are the so called 
4-, 6-, and 3-node reconstruction algorithms, whereby the node indications 

and the coordinate systems of the particular algorithms are given in Fig. 2. 
The correlation between the collected charges Qn at the particular nodes n 
and the reconstructed event positions {u,v) can be realized for the particular 
algorithms in the x- and y-direction as follows: 
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4- node algorithm (events somewhere in cell 1; x^y e [—§,§]) 



1^4 = — With Qxy4. = Qi + Q2 + Q4 + Q5 

(1) 

9 {Q4 + Q5) - {Q1 + Q2) 
Vi — - ■ ^ 

2 Qxy4 

6-node algorithm (node 1 and 4 (node 4 and 5) carry maximum signal for 
reconstruction in x-direction (y-direction) ; x e [0, |] and y e [— f ,0]): 



Uq^ 9 with Q^Q = Qo + <5i + Q2 + Q3 + Q4 + Q5 



Qx6 
)-( 



(2) 



3- node algorithm (node 4 carries maximum signal; x,y e [— f ' 2])• 



^i3 = 5' • with (5a;3 ^ Qs + Q^ + Qb 



Qx3 

7 — 
Qy3 



9- with Qy3 ^ Q1 + Q4 + Q7 



(3) 



The sum of the collected charges in x- and y-direction of the particular al- 
gorithms are abbreviated with Qxy4, Qx6, Qy6, Qx3 and Qy^. The higher the 
ratio of Rl / R2 the more pronounced is the cell like character and the electrical 
screening among the cells. This leads to a treatment of the readout structure 
consisting of strongly separated cells which makes the application of the 4- 
node algorithm stressing this cell structure most favourable. In practise the 
charge screening of the cells is always finite. Therefore, the application of other 
algorithms featuring symmetry planes at the cell border (6-node algorithm) 
and the readout node (3-node algorithm) is sensible. 

The most simplest algorithm imaginable is the 2-node algorithm. But since 

this algorithm only takes two nodes into account it is strongly affected by sys- 
tematic effects (like cross-talk or variations of the amplification of the pream- 
plifiers). Therefore, the 2-node algorithm, when applied in practise, is inferior 
to the other algorithms and is not further considered. 
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Fig. 3. Reconstructed positions within a 8 x 8mm^ cell {Rl = lOOkfi/D and 
R2 = 2kri/n) by means of a) the 4-node algorithm, b) the 6-node algorithm and 
c) the 3-node algorithm. 

3.3 Reconstruction properties of the linear algorithms 
3.3.1 Position reconstruction 

By means of the 4 x 4-cell diffusion simulation it is possible to investigate 
the reconstruction behaviour of the three algorithms. A series of individual 

simulations has been carried out impinging a short driving current successively 
at every grid point (homogeneous grid point spacing Ax = Ay = 200 /im). 
All positions are calculated after the 25 nodes of the 4 x 4-cells have collected 
almost 100 % of the impinged charge. Altogether one obtains therefore 41 x 41 
possible positions within one single cell. 

In Fig. 3 the reconstructed positions of the 4-, 6- and 3-node algorithm are 
shown within a 8 x 8mm^ cell for R1/R2 — 50. The 4-node algorithm shows 
large distortions close to the cell borders which means large deviations between 
the elements of the two spaces {x, y) and {u, v), whereas this algorithm is nearly 
distortion free in the centre of a cell. On the other hand both the 6-nodc 
and the 3-node algorithm have their maximum distortions in the cell centre. 
Close to the cell borders they are more favourable for position reconstruction 
since their distortions are much smaller compared to the 4-node algorithm. 
Obviously, the reconstruction behaviour of the 6- and the 3-node algorithm is 
quite similar, since they use comparable symmetry planes and readout nodes. 

For a homogeneous input illumination - we impinged charges on every grid 
point - a perfect reconstruction algorithm would result in a homogeneous 
output population density. Since none of the introduced algorithms is a per- 
fect reconstruction algorithm (s. Fig. 3) variations in the population density 
become obvious. Mathematically the correlation between the true positions 
{x{u,v),y{u,v)) and the reconstructed positions {u,v) is given by a two- 



7 



dimensional Jakobi-determinant J{u,v) (coordinate transformation), which 
correlates the two population densities B{x,y) and B{u,v): 

dx dy B(x,y) = j dudv J{u,v) ■ B(^x{u,v),y{u,v)^ = jdudvB(u,v) (4) 



Therfore, we obtain 



B{u, v) — J{u, v) ■ b(x{u, v),y{u, v)^ — 



dx{u,v) dy{u,v) 



du du 
dx{u,v) dy{u,v) 



dv 



dv 



B{x,y) (5) 



with the Jakobi-determinant J: 



J 



dx dy 
du dv 



dx dy 
dv du 



du dv 
dx dy 



du dv 
dy dx 



(6) 



whereby the second relationship for J can be derived by the back-trans- 
formation B{x,y) = J{x,y) ■ B{u{x,y),v{x,y)). The derivatives with respect 
to the true positions dx and dy are easier to calculate in a numerical fashion 
since the positions {x, y) are placed on an equidistant lattice in the case of a 
homogeneous illumination. Therefore the second relationship for J from Eq. 
(6) is used for further numerical calculations. 

A homogeneous illumination like in this simulation means B{x,y) — const, 
which is set in the following (arbitrarily) to B(x, y) — 1. This leads to B(u, v) — 
J{u,v) (Eq. (5)). In Fig. 3 the relationship between the true positions {x,y) 
and the reconstructed positions {u, v) is shown for quantised distances. It is 
obvious that J is unequal to one which directly leads to an inhomogeneous 
population density B{u, v). This is also visible by the variation of the density 
of the black points in Fig. 3. A higher density of points {u, v) is equivalent to a 
higher population density B(u, v). In Fig. 4 the population density B{u, v) for 
the 4-nodc algorithm is shown as an example. A perfect algorithm which con- 
nects the reconstructed positions with the true positions x{u, v) and y{u, v) by 
x{u) — u and y{v) ~ v would lead to a population density B{u, v) = B{x, y). 

Unfortunately, it is a non-trivial task to find a suitable transformation between 
the two spaces {x.y) and {u^v) since the functions x{u^v) and y{u,v) are ex- 
plicitly needed. Starting with one of the linear algorithms (e.g. the 4-node al- 
gorithm) one will end up in major problems: Since all the data were calculated 
numerically with a finite grid spacing (e.g. 200 /xm) one has to approximate 
the functions x{u,v) and y{u,v) e.g. by polynomials or Fourier series to get 
an analytical approach. Thereby, it showed up that one has to take too many 
orders into account to decrease the deviation between positions reconstructed 
by the simulation and the polynomials. But even if there exists a satisfying 
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R1/R2=100/2, 4-nodealgor. 




Fig. 4. Population density B{u,v) plotted for the 4-node algorithm for 
Rl = lOOkO/D and R2 = 2kO/n obtained by Eq. (5) with a homogeneous il- 
lumination B{x,y) = 1. 

solution of an analytical expression for x{u, v) and y{u, v) another much more 
serious problem arises. One has to take carefully the spatial resolution into 
account. As it will be seen later on, the transformation for x{u, v) and y{u, v), 
which will correct the reconstructed positions {u,v), will dramatically deteri- 
orate the spatial resolution at the positions where larger corrections have to 
be apphed. This will be for example the case at the cell borders when using 
the 4-node algorithm where we have larger population densities (s. Fig. 4) and 
therefore larger position corrections. 

Also the access via a direct measurement is closed since an extensive mechan- 
ical scan has to be done; e.g. by means of a fine needle impinging a certain 
amount of charge onto the readout structure. At least one has to scan - de- 
pending on the desired accuracy - a few thousand positions which will lead 
to an enormous effort. Besides, the true positions {x, y) are hard to determine 
due to the finite mechanical accuracy of the setup. 

3.3.2 Spatial resolution 

In most position sensitive detector systems the spatial resolution represents 
an important parameter which has to be optimised depending on the desired 
application. Since in interpolating systems like in our case the reconstruction 
algorithms themselves affect the spatial resolution additionally (besides e.g. 
the signal-to- noise ratio (SNR) or the pixel-size for pure pixel detectors), the 
infiuence of the algorithms has to be analysed separately in this respect. 

The absolute values of the spatial resolution of the different algorithms are 
determined by the SNR. Since the spatial resolution of all algorithms is oc 
SNR~^ (s. b.) the absolute noise resistance it! is not exphcitly needed when 
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the various algorithms are compared to each other. In the following we treat 
the parallel noise of the readout structure as the dominant noise source, which 
mainly determines the spatial resolution. Every readout node collecting the 
charge Qi is subjected to Nyquist- Johnson noise 



AQ(,) = JiM^, (7) 



with t denoting the charge integration time. Thereby, the collected charge 
including the noise contribution is normally distributed around Qi{t). For the 
total spatial resolution the electronic noise would also have to be considered. 
It has deliberately not been included in the calculations to leave the discussion 
independent from the electronics used. For most practical cases it will be in 
the same order of magnitude or smaller than the noise contribution, presented 
here. 

Taking into account the particular algorithms (Eqs. (l)-(3)) and assuming a 
constant gaussian charge noise AQ (Eq. (7)) error propagation leads to the 
spatial resolution in ti-direction: 



where (x, y) and t indicate the true event position (impact position) upon the 
readout plane and the charge integration time, respectively. The calculation 
in t;-direction behaves correspondingly. The sum runs over all participating 
nodes of the particular algorithms. Since some charges Qi appear both in the 
u- and the v-coordinate the two coordinates are correlated. Therefore, the 
noise causes a spatial smearing which can be expressed by a two-dimensional 
gaussian distribution with q describing the strength of the correlation (e.g. 
[16]): 



E{x,y) 



1 




X — u) {y — v) 2g(x — u){y — v) 



+ 



(9) 



For the 4-, 6- and 3-node algorithms one obtains for the spatial uncertainty 
(corresponding coordinate systems s. Fig. 2): 
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4-node algorithm: 

au4{u4t),t) ^ ^J4u,{ty + • ^^^^^ = fuA{u4{t),g) ■ (SNR,,4)-' 

COv(m4, V4, t) = 4U4V4 = 4M4t'4 • (SNR3;y4)~^ 

AU4V4 

g{u4, V4, t) 



^{4.ul + g^){4vl + g') 

(10) 



6-node algorithm: 



<Tu6{u6{t),t) ^^6ue{ty + 4g^- ^^^^\ = U{u^{t),g) ■ (SNR^g)"! 
^ Qx6{x,y,t) 

a,,{v,{t),t) = ^/Mtj^+V • = fv6{v,{t),g) ■ (SNR,6)-' 

(Ag)2 (111 

cov(m6, ^6, = (4m6V6 + 25^(^6 - ve) - g 



g{u6,VG,t) 



AuqVq + 2g{uQ -Vq) - g 



2 



^J{6ul + Ag^){6vl + Ag^) 
3-node algorithm: 

C7u3iu3it),t) ^^3us{ty + 2g^- ^^^^\ = fu3{us(t),g) ■ (SNR^a)"' 

^!x3\x, y, t) 

ay3{v3{t),t) - ^J3vs{ty + 2g'- = fv3{v3{t),g) ■ (SNR^a)"^ 

COy{U3,V3,t) = U2V3 



(12) 



Qx3 Qy3 

. . Mat's 



^(3«i + 2^2)(3^2 + 2^2) 

The cell size is denoted by Again the reconstructed positions are indicated 
by (u(x,y),v(x,y)) together with a number as a subscript for the particular 
algorithm. The terms of the spatial resolutions can be separated into a geomet- 
rical factor / and a factor consisting of the SNR. Fig. 5 shows schematically 
the results of the spatial resolutions du — 2.355 • (T„ (fwhm) as a function of 
the distorted reconstructed position («„, v„) for all three algorithms within 
a single 8 x 8mm^-cell taking both the geometrical factor / and the SNR- 
factor into account. Obviously, the 4-node algorithm shows up to be the most 
favourable algorithm since it provides the best spatial resolution. Since the 
3-node algorithm collects in the cell centre only about half of the impinged 
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R1/R2=100/2 




Fig. 5. Spatial resolution du = 2.355 • (T„ (fwhm) for the 4-, 6- and 3-node al- 
gorithm as a function of the reconstructed positions {un,Vn) for Rl = lOOAjfi/D 
and R2 = 2kQ/0 after the complete charge has been integrated (t = 500ns). The 
driving current contains a charge equivalent of about 1.7- 10^ e~. The parallel noise 
resistance of the readout structure has been numerically calculated to i? Ri 13.6 kJ7. 

charge the SNR-factor decreases and therefore the spatial resolution deteri- 
orates in comparison to the other algorithms. Only at the nodes the spatial 
resolution for the 3- and the 4-node algorithm becomes similar, since the geo- 
metric factors / = V2g (cf. Eqs. (10) and (12)) and the SNR-factors are the 
same. 



Actually, the spatial resolution in Fig. 5 is plotted as a function of the recon- 
structed positions {u,v) and not of the true positions {x,y). If one corrects 
the reconstructed positions by means of a coordinate transformation which 
transforms {u, v) {x, y) the spatial resolution of the corrected positions Ux 
can be estimated by: 



y)}^ ^ 



dudv du dv\'^ 



dxdy dydxj 

Eq. (13) determines the spatial resolution of the corrected positions as a 
function of the covariance, the spatial resolutions and of the applied 
algorithm (cf. Eqs. (10)-(12)) and the derivatives between the coordinates 
(m, w) and {x,y). In this equation all appearing {u,v) should be treated as 
{u{x,y),v{x,y)). Fig. 6 shows exemplarily the numerically calculated spatial 
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Fig. 6. Spatial resolution dx = 2.355 ■ ctx (fwhm) for the 4-node algorithm as a 
function of the true positions {x,y) for Rl = W0kQ,/O and R2 = 2kQ,/0. The 
graph is plotted for the same parameters like in Fig. 5. Since both the 6- and the 

3- node algorithm show quite the same behaviour close to the cell borders only the 

4- node algorithm is plotted exemplarily. 

resolution dx = 2.355 • (fwhm) obtained by Eq. (13) as a function of the 
true positions (x, y) for the 4-node algorithm. Since at the cell borders where 
the population density of all algorithms is relatively high (s. Fig. 3) larger 
corrections have to be applied the spatial resolution deteriorates by almost a 
factor of 3-4 for all algorithms. 



3.4 Introduction of the optimised 463-node algorithm 

In the following we present a possible solution which takes care of both the 
correct position reconstruction and the spatial resolution, by combining the 
4-, 6- and 3-node algorithms to one 463-node algorithm. 

All algorithms converting a homogeneous population density B{x,y) — 1 into 
an inhomogeneous population density B{u,v) 7^ 1 which is equivalent to a 
wrong position reconstruction are not favourable. Even for the case of a later 
position correction these algorithms will always suffer from a bad spatial reso- 
lution at the positions where larger position corrections have to be applied (cf. 
end of Sec. 3.3.2). Therefore, it would be much more advantageous to find an 
algorithm which is capable to reconstruct almost the correct positions which 
require only small subsequent position corrections and which therefore only 
weakly influence the spatial resolution. The idea of combining the 4-, 6- and 
3-node algorithm to the 463-node algorithm shows up to be a possible solution 
fulfilling almost the requirements mentioned above. 
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Fig. 7. Choice of the mixing matrices Ax and B^. The first bifurcation mixes 
between the 4-node and either the 6- or the 3-node algorithm. The criterion for 

aij G [0, 1] is an optimised position reconstruction which will be described in the 
text. The second bifurcation distinguishes between the 6- and the 3- node algorithm. 
The criterion for bij G {0, 1} is given by the spatial resolution au6 and a^s- 

Based on the results of the simulation model giving a relationship between 
{u, v) and {x, y) the combination of the algorithms in a;-direction is done by 
mixing matrices — ((aij)^) and = {{bij)^), defined for every grid point 
in the 8 x 8mm^ cell. The elements a^j e [0,1] mix the 4-node algorithm 
continuously with either the 6-node or the 3-node algorithm, whereby the 
elements of the second mixing matrix bij G {0, 1} consisting only of the values 
or 1 choose either the 6-node {bij = 1) or the 3-node {bij = 0) algorithm. The 
following equation defines the complete 463-node algorithm, whereby Ux — Oij 
and bx — hj are those elements of the matrices corresponding to the true 
position {x,y). The origin of the 4-node algorithm is now shifted to readout 
node 4 (cf. Fig. 2): 



For symmetry reasons the relationship between the mixing matrices in x- 
and y-direction is simply: Ay = A^ and By = B^. The procedure for the 
combination of Ax and Bx is shown in the schematic diagram in Fig. 7. In 
the following the derivation of the matrices Ax and Bx is discussed for the 
x-direction. The y-direction behaves correspondingly. 

The spatial resolution a^e and a^s of the 6- and the 3- node algorithm is 
the criterion for the matrix Bx with bij G {0,1} (cf. Eqs. (11) and (12)). If 
Cu6 < CTus the 6-node algorithm is used {bij = 1); if au6 > o'us the 3-node 




, /, (Q2 + Q5)-(go + Q3) 

{^-ax)9- <bx ^ 



+ (1 - bx) ■ 



Qx3 



] 

(14) 
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algorithm is used {bij = 0). Since the deviations between the reconstructed 
positions of the 6-node and the 3-node algorithm are small (s. Fig. 3) also the 
influence of the points of discontinuity {bij = bij = 1) is expected to be 
negligible. 

The criterion for the combination between the 4-node and the 6- or the 3- node 
algorithm is an optimised position reconstruction which means a minimum 
deviation A between u and x. Together with ^6,3 = bxUQ + (1 — bx)us we get 
as a criterion: 



min A = min 



X - 1i463 



= min I (a: - axU4 - (1 - ax)ue,3) | , 



(15) 



whereby again = and b^ = bij. The minimum condition leads to 
dA/dax = 0. Since d'^A/da^ > the criterion for a minimum is fulflUed. 
As a result we obtain for the matrix elements Ux' 

X — uqs I if > 1 then ax is set to 1 , 

ttx — '— < (16) 

1*4 — 1*6,3 I if < then ax is set to 

The constraint 

< < 1 , (17) 

has to be set because for some positions mixing factors \ax\ ^ 1 would end 
up in additional distortions due to noise or systematic effects (cf. Sec. 4). 
Consequently, this constraint leads to the fact that some positions {u,v) are 
still not reconstructed correctly. 

A graphical illustration of the matrices Ax and Bx arc shown in Figs. 8 and 9. 
The 4-node algorithm determines already the optimum positions in x-direction 
close to the low resistivity cell borders (except close to the readout nodes), 
therefore it is the only algorithm applied (ax = 1). In y-direction close to the 
cell borders the 4-node algorithm becomes unfavourable (cf. Fig. 3 a)) and 
the other participating algorithms are more emphasised (a^: ~ 0). As expected 
the 6-node algorithm offers the best spatial resolution close to the symmetry 
axis of the cell, whereas the 3-node algorithm becomes more favourable in 
x-direction close to the low resistivity cell borders and the readout nodes 
(cf. also Fig. 5). Fig. 10 shows the reconstructed positions obtained with the 
463-node algorithm by means of Eq. (14). The slight distortions for some 
positions mainly close to the nodes are due to the constraint in Eq. (17) 
and can be directly revealed in the population density plot (Fig. 11). A very 
slight inhomogeneity is also visible close to the low resistivity strips at the cell 
borders. In comparison with the population density of the 4-nodc algorithm 
(Fig. 4) the population density of the 463-node algorithm is nearly flat and 
especially close to the cell borders and the readout nodes improved by a factor 
2-4. 
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Fig. 8. Matrix with a.^ £ [0, 1] calculated with Eq. (16) for Rl = lOOfcQ/D and 
R2 = 2 kQ/O [uij = 1 corresponds to the full usage of the 4-node algorithm, aij = 
corresponds to the pure use of either the 6- or the 3-node algorithm) . 
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Fig. 9. Matrix with G {0, 1} plotted for Rl = 100A;J7/n and R2 = 2k9./U 
{bij = 1 (black area) corresponds to the full application of the 6-node algorithm and 



(white area) corresponds to the full application of the 3-node algorithm). 



Besides this as a consequence also the spatial resolution is improved. In Fig. 
12 the spatial resolution of the 463-node algorithm is plotted as a function of 
the reconstructed positions {u,v). In fact the error propagation (Eq. (8) and 
Eq. (14)) leads for the 463-node algorithm to a non-trivial equation for the 
spatial resolution. However, as expected due to the mixing matrices and 
Bx close to the nodes and the cell borders in v-direction the spatial resolution 
is comparable to that obtained by the single 3-node and 6-node algorithm (cf. 
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Fig. 10. Reconstructed positions by means of the 463-node algorithm within a 
8 X 8mm2 cell {Rl = lOOkJ^/D and R2 = 2kJ2/n). The corresponding mixing 
matrices and are shown in Figs. 8 and 9, respectively. 
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Fig. 11. Population density B{u,v) plotted for the 463-node algorithm for 
Rl = lOOkJl/D and R2 = 2kO/n obtained by Eq. (5) with a homogeneous il- 
lumination B{x,y) = 1. 



Fig. 5). Close to the cell borders in it-direction the behaviour of the single 4- 
node algorithm is reflected. Since the positions are indeed well reconstructed 
(Fig. 10) the spatial resolution of the 463-node algorithm as a function of the 
true positions (x, y) is very similar to that of the reconstructed positions {u, v) 
(Fig. 12). Only close to the nodes the spatial resolution is slightly deteriorated 
caused by the small position corrections (s. Eq. (13)). 
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Fig. 12. Spatial resolution du = 2.355 • o"„ (fwhm) for the 463-node algorithm as a 
function of the reconstructed positions {u,v) for Rl = lOOkQ/O and R2 = 2kQ/0. 
The graph is plotted for the same parameters like in Fig. 5. 

3.5 Application of the 463-node algorithm on measured data 

The application of the 463-node algorithm (Eq. 14) demands for the mixing 
matrices A and B. Since these matrices are actually a function of the true po- 
sitions (x, y) which are not known a priori in a measurement we have to apply 
an iterative technique. To obtain a first approximate position an assumption 
for suitable elements a^j and bij is needed. For that reason the reconstructed 
position is calculated firstly both by means of the 6-node and the 4-node al- 
gorithm. At a distance of c? = 0.1 g (with g as the cell size) around the cell 
borders the 6-node algorithm (uq,vq) is used for a first estimation of i,j to 
obtain the matrix elements a^j and 6^. Else the 4-node algorithm {u4,V4) is 
used for the estimation of and bij. Then, the 463-node algorithm posi- 
tion is calculated for the first time, which gives new positions and therefore 
an estimation of new aij and bij. This iterative procedure is repeated several 
times until the reconstructed positions (^403,^403) converge. Since the diffu- 
sion simulation gives matrices A and B with 41 x 41 elements corresponding 
to a grid spacing of 200 /^m (Figs. 8 and 9) we have used an interpolation 
routine of IDL [17] to achieve a finer resolution. If the virtual pixel size is e.g. 
chosen to 200 x 200 fira'^ we use finer matrices A and B with 81 x 81 elements, 
corresponding to 100 fim grid spacing to allow a better convergence of the 463- 
node algorithm. Since this method converges at almost all positions within the 
cell this iterative loop is used for all image reconstructions of measured data 
presented in the following. 

Fig. 13 shows the response of the inner 5x5 cells (= 40 x 40 mm^) of a 
detector with a PCB-readout structure to an uniform illumination {B{x,y) ~ 
1) reconstructed with the 4-node algorithm (Eq. 1). A ^^Fe-source has been 
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Fig. 13. Reconstructed image of a flatfield illumination using the 4-node algorithm. 
A virtual pixel size of 200 x 200 ^vo? is chosen. The image contains about 4.57 • 10^ 
photons corresponding to a mean number of photons per pixel ~ 114. The stan- 
dard deviation of the intensity per pixel distribution amounts to 54.9. 



used for illumination {E^ = 5.9 keV). The detector has been operated with an 
Ar/C02 (90/10) gas filling at standard pressure. 



The low resistivity cell borders of the PCB-structure have a width of about 
(170 ± 5/im). The special resistive materials used for the silk-screen print- 
ing process have surface resistances of Rl — lOOkQ/D and R2 — IkQ/D, 
respectively. Unfortunately the printing and burning processes have a strong 
influence on the resistances and the exact ratio of -Rl / R2 is hardly predictable. 
Since the depletions at the cell borders when using the 4-node algorithm show 
a good regularity in all cells (Fig. 13) we conclude that the relative accu- 
racy of the surface resistances is in the percent range. The resulting image 
when using the 463-node algorithm is shown in Fig. 14. The matrices A and 
B are created for Rl = lOOkfi/D and R2 = lOkfi/D since this ratio sim- 
ply fits best, whereby the distortions of the 4-node algorithm depend less 
than linear on the ratio R1/R2. This means that one obtains e.g. only lit- 
tle changes concerning the distortions when treating a ratio of R1/R2 — 10 
instead of R1/R2 — 20. The population density -B(m463, '^463) becomes much 
more homogeneous in comparison to the 4-node algorithm. In addition the 
standard deviation of the intensity per pixel is dramatically reduced from 
(Jiv = 54.9 to (Tiv = 16.2 which is indeed much closer to the theoretical limit 

of (Tiv-Poisson = V^ni 10.7. 
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Fig. 14. Reconstructed image from the same data as in Fig. 13 using the 463-node 
algorithm. A virtual pixel size of 200 x 200 iivo? is chosen. The standard deviation 
of the intensity per pixel distribution amounts to 16.2. The dotted circles indicate 
systematic effects. Uppermost circle: not properly working preamplifier. Middle cir- 
cle: low resistivity strip is interrupted. Lowermost circle: defect in a copper layer in 
one of the three GEM structures used for gas amplification. 

4 Non-linecir corrections 

Although the reconstructed image with the 463-node algorithm (Fig. 14) rep- 
resents a mentionable progress compared to images reconstructed e.g. with the 
4-node algorithm (Fig. 13) still some inhomogeneities in the population den- 
sity are obvious. These artifacts can mainly be attributed to systematic effects 
(which partly occur during the measurement). As already described, even in 
the theoretical case the 463-node algorithm leads to a slight overpopulation 
close to the cell borders and especially around the nodes (cf. Fig. 11). System- 
atic electronic effects, like gain variations of the preamplifiers and cross-talk in 
the preamplifiers/cables can possibly influence the image. Besides this, a too 
short integration time of the preamplifiers can lead to an overpopulation (as 
can be shown by simulations) of the cell borders. Further systematic effects 
stem from the readout structure itself, like the finite dimensions of the readout 
nodes, the finite accuracy of the width of the low resistivity strips and local 
inhomogeneities of the surface resistances. 

4-1 Homogeneous population density 

In order to decrease the remaining inhomogeneities around the readout nodes 
and the cell borders (cf. Fig. 14) an ansatz is made for x{u,v) and y{u,v), 
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Fig. 15. Measured averaged population density B{u, v) of the 5x5 cells of Fig. 14. 
The average value of the population density is normalised to unity. 

whereby the abbreviation u — u^^^ and v — V4q^ is used in the following. 
Since we know that the incident population density B{x, y) for a homogeneous 
illumination is equal to unity Eq. (5) reads as follows: 

B(x.y) = l= . /'"'^^ . ^ B(„,„) = |£|^ (18) 

^ dx_dy_ _d^dy_ ^ ' ' dudv dv du 

du dv dv du 

The functions ,x'(?i. v) and y{u,v) have now to be optimised in such a fashion 
that the derivatives duX d^y — d^x duy correspond to the measured population 
density B{u,v). Fig. 15 shows the measured mean population density B{u,v) 
obtained by the superposition of the 5x5 cells shown in Fig. 14. The mean 
value of B{u,v) is normalised to unity; also the cell size g is normalised to 
one {u,v G [0,1]). We have composed the coordinate transformation func- 
tions x{u,v) and y{u,v) by suitable combinations of different one- and two- 
dimensional Gaussian functions with the cell centre as the point of symmetry. 
The final coordinate transformation functions x{u,v) and y{u,v) have alto- 
gether 6 independent parameters. Fig. 16 shows the population density as a 
result of the minimisation of \B{u. v) — [duxd^y — dyxduy)\'^- The measured 
population density B{u,v) (Fig. 15) and the population density obtained by 
the coordinate transformation functions x{u,v) and y{u,v) (Fig. 16) are in 
good agreement. The knowledge of x{u,v) and y{u,v) allows now to correct 
for the small inhomogeneities around the nodes and the cell borders. By this 
procedure, the reconstructed positions of the 463-node algorithm {u, v) are 
transformed into the corrected positions x{u,v) and y{u,v). Fig. 17 shows the 
resulting non-linear corrected image. The standard deviation of the intensity 
per pixel is further reduced from = 16.2 (Fig. 14) to = 13.6 which is 
even closer to the theoretical limit of c at— Poisson — 

^/lU^ 10.7. 



21 




Fig. 17. Reconstructed image of the flatfield illumination using the 463-node algo- 
rithm with subsequent non- linear corrections. The standard deviation of the inten- 
sity per pixel distribution amounts to 13.6. 

4-2 Inhomogeneous population density 

Since all assumptions made in the previous sections are based on a homoge- 
neous population density B{x,y) = 1 we have to test the 463-node algorithm 
(Sec. 3.4) and the corresponding non-linear corrections separately with an 
inhomogeneous population density. This can e.g. be performed with a suit- 
able collimator. More images (even time-resolved measurements) in the field 
of biological and chemical diffraction samples recorded with the ViP-detector 
prototype can be found in Ref. [12]. 
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Fig. 18. On the left hand side the image of the "SAXS" aperture recorded with a 
PCB-readout structure is shown. On the right hand side a photo of the aperture 
taken with a standard digital camera against the sunlight is depicted. The distance 
between detector and source amounted to about 2.5 m, therefore parallax is very 
small. The illuminated spots at the bottom correspond to holes in the aperture with 
diameters of 280, 380, 480, 580 and 680 fim. The aperture is slightly tilted by an 
angle of about 0.6°, which can be recognised for example by the vertical pixel jump 
(pixel quantisation) at the bottom of the "2" (dashed-dotted line) . The sizes of both 
depicted images amount to 4.4 x 4.4 cm^. 

Fig. 18 shows the image of a laser cut 1mm thick stainless steel aperture 
containing "SAXS 2D" letters (Small Angle X-ray Scattering) and five holes 
with increasing hole diameters. The image has been recorded with a PCB- 
rcadout structure using photons of an energy of 6.4 keV (fluorescence of a Fe- 
target in a 8keV synchrotron beam). The detector was filled with a Xe/C02 
(90/10) gas mixture at a pressure of 1.3 bar. The image has been reconstructed 
with the linear 463-node algorithm with subsequent non-linear corrections 
at the nodes and at the cell borders. We performed an additional flatfleld 
correction with a flatfleld image recorded in our laboratory one week after the 
aperture measurement with a different gas flUing and a different X-ray source. 
Since this flatfleld image is different from one which we would have obtained at 
the beamline directly after the aperture measurement (with identical detector 
parameters), the corrected image still shows up some artifacts like visible areas 
around the readout nodes (marked with dotted lines) which should have been 
suppressed by the correction of the proper flatfleld image. 

However, the recorded image compares well to the image of the aperture (right 
hand side image in Fig. 18). Only the middle part of the second "S" (dashed 
line) looks slightly distorted; at this readout channel the preamplifler was not 
working optimally (compare to uppermost dotted circle in Fig. 14). Also the 
reproduction of the holes shows a good agreement. 
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5 Conclusion 



The simplest linear algorithms do not reconstruct the event positions correctly. 
This leads to distortions since the true positions {x, y) and the reconstructed 
positions {u,v) differ up to a certain extent. When these algorithms are cor- 
rected afterwards by a suitable coordinate transformation connecting the two 
spaces {u, v) and {x, y) the spatial resolution becomes worse especially next 
to the cell borders and the readout nodes. 

Therefore, the 463-node algorithm is introduced. This complex algorithm con- 
sisting of combinations of linear algorithms is a possible solution for posi- 
tion reconstructions using the two-dimensional interpolating resistive readout 
structure described above. This algorithm has been optimised with respect to 
both the correct position reconstruction and the spatial resolution. Indeed, 
the reconstruction of a homogeneous population density as a measure for the 
quality of an image has been considerably improved - even close to the Poisson 
limit - when the 463-node algorithm with subsequent non-linear corrections 
was applied. Nevertheless, some inhomogeneities mainly caused by systematic 
effects arc still visible in the images. Since these disturbing effects can not be 
corrected afterwards they should be kept as small as possible. 



References 

[1] A. Bressan, R. De Olivcira, A. Gandi, J.-C. Labbe, L. Ropelewski, F. Sauli, 

D. Mormann, T. Miillcr, H. J. Simonis, Two-dimensional readout of GEM 
detectors, Nucl. Instr. and Meth. A425 (1999) 254 261. 

[2] S. Baehmann, S. Kappler, B. Ketzer, T. Miiller, L. Ropelewski, F. Sauli, 

E. Schulte, High rate X-ray imaging using multi-GEM detectors with a novel 
readout design, Nucl. Instr. and Meth. A478 (2002) 104-108. 

[3] J. P. Cussonneau, M. Labalmc, P. Lautridou, L. Luquin, V. Metivier, 
A. Rahmani, T. Reposeur, 2D localization using resistive strips associated to 
the Micromegas structure, Nucl. Instr. and Meth. A492 (2002) 26-34. 

[4] R. Lewis, Multiwire Gas Proportional Counters: Decrepit Antiques or Classic 
Performers?, J. Synchrotron Rad. 1 (1) (1994) 43-53. 

[5] H. J. Besch, M. Junk, W. Meifiner, A. Sarvestani, R. Stichlcr, A. H. Walenta, 
An interpolating 2D pixel readout structure for synchrotron X-ray diffraction 
in protein crystallography, Nucl. Instr. and Meth. A392 (1997) 244-248. 

[6] F. Sauli, GEM: A new concept for electron amplification in gas detectors, Nucl. 
Instr. and Meth. A386 (1997) 531-534. 



24 



[7] Y. Giomataris, P. Rebourgeard, J. P. Robert, G. Charpak, MICROMEGAS: 
a high-granularity position-sensitive gaseous detector for high particle-flux 
environments, Nucl. Instr. and Mcth. A376 (1996) 29-35. 

[8] F. Bartol, M. Bordessoule, G. Chaplier, M. Lemonnier, S. Megtert, The C.A.T. 
Pixel Proportional Gas Counter Detector, J. Phys. Ill France 6 (1996) 337-347. 

[9] A. Sarvcstani, H. J. Besch, M. Junk, W. Meifiner, N. Sauer, R. Stiehler, A. H. 
Walenta, R. H. Menk, Study and application of hole structures as gas gain 
devices for two dimensional high rate X-ray detectors, Nucl. Instr. and Meth. 
A410 (1998) 238-258. 

[10] A. Orthen, H. Wagner, H. J. Besch, R. H. Menk, A. H. Walenta, 
U. Wcrthenbach, Investigation of the performance of an optimised MicroCAT, 
a GEM and their combination by simulations and current measurements, Nucl. 
Instr. and Meth. A500 (2003) 163-177. 

[11] A. Orthen, H. Wagner, H. J. Besch, S. Martoiu, R. H. Menk, A. H. Walenta, 
U. Werthenbach, Gas gain and signal length measurements with a triple-GEM 
at different pressures of Ar-, Kr- and Xe-based gas mixtures, Nucl. Instr. and 
Meth. A512 (2003) 476-487. 

[12] A. Orthen, H. Wagner, S. Martoiu, H. Amenitsch, S. Bernstorff, H. J. Besch, 
R. H. Menk, K. Nurdan, M. Rappolt, A. H. Walenta, U. Werthenbach, 
Development of a two-dimensional virtual pixel X-ray imaging detector for time- 
resolved structure research, submitted to J. Synchrotron Rad. (2003). 

[13] A. Sarvestani, N. Sauer, C. Strietzel, H. J. Besch, A. Orthen, N. Pavel, A. H. 
Walenta, R. H. Menk, Microsecond time-resolved 2D X-ray imaging, Nucl. Instr. 
and Meth. A465 (2001) 354-364. 

[14] S. Martoiu, A. Orthen, H. Wagner, H. J. Besch, R. H. Menk, K. Nurdan, A. H. 
Walenta, U. Werthenbach, Intelligent local trigger technique for a multi-cell 2D 
interpolating resistive readout, to be submitted to Nucl. Instr. and Meth. A 
(2003). 

[15] H. Wagner, H. J. Besch, R. H. Menk, A. Orthen, A. Sarvestani, A. H. 
Walenta, H. Walliser, On the dynamic two-dimensional charge diffusion of 
the interpolating readout structure employed in the MicroCAT detector, Nucl. 
Instr. and Meth. A482 (2002) 334-346. 

[16] S. Brandt, Data Analysis, Spinger, New York, 1999. 

[17] IDL 5.3, Research Systems Inc., Boulder, CO, USA. 



25 



